UNCLASSIFIED 


AD  NUMBER 


AD805554 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  DEC  1966. 
Other  requests  shall  be  referred  to 
Ballistic  Systems  Div.,  AFSC,  Norton  AFB, 
CA. 


AUTHORITY 


SAMSO  USAF  ltr,  28  Feb  1972 


THIS  PAGE  IS  UNCLASSIFIED 


A  SUMMARY  ANALYSIS  OF  LABORATORY  HYPERSONIC 
WAKE  TRANSITION  EXPERIMENTS 


A.  Goldburg 


RESEARCH  NOTE  672 
December  1966 


supported  jointly  by 


BALLISTIC  SYSTEMS  DIVISION 
DEPUTY  FOR  BALLISTIC  MISSILE  RE-ENTRY  SYSTEMS 
AIR  FORCE  SYSTEMS  COMMAND 
Norton  Air  Force  Bote,  California 
under  Contract  AF  04(694) -690 


ADVANCED  RESEARCH  PROJECTS  AGENCY 
monitored  by  tbe 
ARMY  MISSILE  COMMAND 
UNITED  STATES  ARMY 
Redstone  Arsenal,  Alabama 
ondei  Contract  DA -01  -021  -  AMC-  1200ff  (Z) 
(part  of  Project  DEFENDER)  7 


Best  Available  Copy 


HSD-TR-bb-475 


RESEARCH  NOT 


: - A  S liMMARY  ANALYSIS  Or  LAhORATOk  if  HY PElCSONIC 

WAKE  TRANSITION  EXPERIMENTS 

by 

A.  Goldburg 


AVCO  EVERETT  RESEARCH  LABORATORY 
a  division  of 
AVCO  CORPORATION 
Everett,  Massachusetts 


December  1966 


supported  jointly  by 


BALLISTIC  SYSTEMS  DIVISION 
DEPUTY  FOR  BALLISTIC  MISSILE  RE-ENTRY  SYSTEMS 
AIR  FORCE  SYSTEMS  COMMAND 
Norton  Air  Force  Base,  California 

under  Contract  AF  04(694)-690 


ADVANCED  RESEARCH  PROJECTS  AGENCY 
monitored  by  the 
ARMY  MISSILE  COMMAND 
UNITED  STATES  ARMY 
Redstone  Arsenal,  Alabama 

under  Contract  DA-01 -021 -AMC-120G5  (Z) 
(part  of  Project  DEFENDER) 


FOREWORD 


This  report  has  been  supported  jointly  by  Avco  Everett  Research 
Laboratory  under  Contract  AF  04(694)- 690  for  Ballistic  Systems  Division, 
Deputy  for  Ballistic  Missile  Re-entry  Systems,  Air  Force  Systems  Com¬ 
mand,  Norton  Air  Force  Base,  California,  and  Advanced  Research  Projects 
Agency,  monitored  by  the  Army  Missile  Command,  United  States  Army, 
Redstone  Arsenal,  Alabama  under  Contract  DA-01-021-AMC-12005  (Z) 

(part  of  Project  DEFENDER).  The  secondary  report  number  as  assigned 
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ABSTRACT* 


Hypersonic  Wake  Transition  data  over  a  range  of  Reynolds  numbers 
have  been  obtained  over  the- past  several  years  at  many  laboratories..  .la.  the 
present  paper,  the  available  experimental  data  is  structured  into  three  laws: 
The  Law  of  the  Far  Wake,  The  Law  of  the  Near  Wake,  and  The  Caw  of  the 
Interpolation  Regime. — These  three  laws  are  then  represented  on— a  single 
Hypersonic  Wake  Transition  Map. 

The  Reynolds  number  is  selected  as  the  correlating  parameter, 
where  local  conditions,  characterizing  the  state  of  the  gas  at  transition 
are  used  for  the  property  values.  The  correlating  constant  is  determined 
from  the  ballistic  range  data.  For  high  Mach  number  flows,  binary  scaling 
allows  the  representation  of  the  Reynolds  number  as  a  pressure  times  a 
length. 


In  the  Far  Wake,  for  the  first  onset  of  wake  unsteadiness,  the  wake 
velocity  profiles  have  reached  their  asymptotic  behavior,  density  gradients 
have  become  small,  and  details  of  body  shape  have  been  lost.  For  this 
region,  the  characteristic  length  is  represented  as  -n/'CqA  based  on  the  total 
drag  of  the  body.  The  Far  Wake  transition  rule  becomes  */C j)A  =  constant. 

In  the  Near  Wake,  local  hypersonic  phenomena  depending  on  body 
shape  dominate  the  flow  field.  Here  the  characteristic  length  is  the  distance 
from  the  wake  origin  to  the  onset  of  unsteadiness,  As  in  boundary 

layer  flow  this  correlation  shows  a  strong  dependence  on  an  appropriate 
local  Mach  number,  Me  ,  at  the  point  of  wake  transition:  p  •  x^  -  constant, 
a  function  of  M_. 

C 

In  the  Interpolation  Regime,  one  traverses  from  the  Near  to  the  Far 
Wake.  The  strong  favorable  pressure  gradient,  the  strong  density  gradients 
and  the  effect  of  body  shape  becomes  less  important  as  one  proceeds  down¬ 
stream  along  the  Wake  in  the  Interpolation  regime.  Thus,  ideas  derived 
from  the  Navier -Stokes  equation  are  useful  and  lead  with  the  aid  of  the 
ballistic  range  data  to  a  sufficiently  accurate  general  interpolation  curve 
to  connect  the  Far  Wake  Law  with  the  Near  Wake  Law.  The  Interpolation 
curve,  like  the  pw  •  V CpA  Far  Wake  and  the  pw  •  x^  Near  Wake  correla¬ 
tions  is  good  to  a  factor  of  Z  in  the  pressure. 

The  Transition  Map  provides  a  method  for  making  the  transition 
curve,  xtr  as  function  of  pressure  ratio  (altitude),  for  a  given  body  at  a 
given  flight  velocity:  Connect  the  Far  Wake  onset  altitude  for  the  n/CqA 
of  the  body  with  the  Near  Wake  behavior  for  the  correspondingly^ 

of  the  body.  Thus  a  complete  qualitative  Wake  Transition  Map  is  presented 
fo?  a  large  range  of  Mach  numbers  and  body  sizes  and  body  shapes. 


♦This  paper  is  a  summary  of  the  author's  work  on  this  topic  over  the  i 
years  he  was  associated  with  Avco  Everett  Research  Laboratory,  It  i 
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1.  INTRODUCTION* 

The  latent  of  the  present  paper  is  to  provide  a  rational  framework 
for  hypersonic  wake  transition  which  lias  a  solid  foundation  in  fluid 
mechanical  ccncoptn  and  understanding,  into  which  valid  transition  in-  ' 
formation  can  he  logically  placed,  and  which  will  be  usefully  structured 

— for  the  -teohnolcgicai  seeds  which  motivate  the  problem,  —  - -  — 

The  wakes  behind  hyper  velocity  bodies  daring  atmospheric  reentry  • 
and  ealiistic  range  flight  can  be  observed  by  means  of  active  {radar}  or 
passive  { self-luminous )  electromagnetic  radiation.  The  patterns  of  ob¬ 
servables  in  the  hypersonic  wake  may  be  thought  of  as  being  produced 
in  the  following  way:  Physical  chemistry  determines  the  observable 
tracers  and  fluid  mechanics  determines  the  mixing  patterns.  These  two  - 
effects  are  formally  related  through  the  coupling  between  the  fluid  mechani¬ 
cal  equations  of  motion  (continuity  and  momentum);  the  equations  related  to 
the  constituents  of  the  fluid  (species  conservation);  and  the  equations  related, 
to  the  constituents'  quantum  mechanical  properties  (energy);  and  the  equation 
of  state.  The  radiating  and  scattering  properties  of  the  gas  are  strongly  in¬ 
fluenced  by  the  fluid  mechanical  mixing  process  which  in  turn  is  strongly  af¬ 
fected  by  whether  the  mixing  is  being  accomplished  solely  through  molecular 
diffusion  or  whether,  in  addition,  there  is  macroscopic  mass  interchange, 
hi  order  to  recover  the  body  description  Vom  the  wake,  information 
available  in  the  wake  must  be  extracted  and  used.  The  wake  is  a  highly  dif¬ 
fusive  region.  In  diffusion  processes,  information  is  lost  as  time  passes. 

In  contradistinction  to  the  wave  equation,  time  in  the  diffusion  equation  pro¬ 
ceeds  in  only  one  direction.  Hence,  the  attempt  to  reconstruct  the  details 

of  body  mass  and  shape  may  require  attention  to  the  local  distributions  of - 

wake  properties. 

The  requirement  to  describe  local  property  distributions  in  fluid 
mechanical  flow  fields  introduces  new  questions  for  ae rodynami cists.  In 
the  past,  aerodynamics  has  been  concerned  with  predicting  time-aver¬ 
aged  and  space -averaged  forces  and  heat  transfer  at  solid  surfaces. 

<c 

This  report  represents  a  summary  of  the  author's  work  in  the  area  of 
hypersonic  wake  transition  over  the  several  years  he  was  associated  with  the 
Avco  Everett  Research  Laboratory.  Publication  of  this  report  does  not  consti 
tute  Avco  Everett  Research  Laboratory's  approval  of  the  findings  and  con¬ 
clusions  therein.  It  is  being  published  only  for  the  exchange  and  stimulation 
of  ideas. 


Fortunately,  one  could  be  quite  wrong  about  the  local  property  distribu¬ 
tions  and  yet  obtain  very  useful  answers  for  surface  viscous  effects. 

.Now,  however,  .the  requirement  is  for  a  description,  local  in  space  and 
local  in  time,  of  property  distributions  of  interest  -  ,  e.  g.  ,  a  velocity  or  a 
particular  species  concentration  —  as  a  function  of  reentry  conditions. 
Further,  if  the  flow  is  time  dependent,  a  description  of  the  local  spec¬ 
tral  distribution  of  the  fluctuations  of  each  important  wake  property 
may  be  required.  In  the  case  of  radar  scattering  by  the  ionized  turbu¬ 
lent  hypersonic  flow  field,  acting  as  a  collection  of  single  scatterers, 
the  spatial  correlation  function  of  the  refractive  index  fluctuations  is 
required.  Further,  the  scattering  is  determined  by  a  narrow  portion 
of  the  spectrum  around  f  k^-kj^J  ),  where  k  is  the  wave  number  and 
X.  the  vector  from  the  illuminating  source;*  and  hence,  its  prediction  is 
sensitive  to  the  accuracy  with  which  such  spectra)  distributions  are  de¬ 
scribed. 

Past  aerodynamic  transition  problems  have  been  concerned  with 
body  related  ohear  layers  -  either  along  or  immediately  separated  from 
the  surface  boundary.  In  such  body-related  shear  layers,  the  hypersonic 
transition  process  from  laminar  to  small  scale  turbulence  is  com¬ 
pressed  into  a  small  distance  relative  to  surface  length  of  interest  (say 
an  airfoil  chord  length)  due  to  the  continuous  generation  of  vortical 
energy  along  the  surface.  The  past  technological  requirement  for  only 
space  end  tjiT=e  raged  surface  information  and  the  narrowness  of  the 
transition  region  relative  to  the  surface  extent  led  to  the  concept  of, 

...  . it  was  *'x  '*  that  the  forces  or  heat  transfer  went  from  the 

transition  tr 

laminar  value  to  the  more  or  less  turbulent  value. 

Transition  in  the  hypersonic  wake  is  more  gradual  since  the  vortical 
energy  is  created  near  the  body  alone  and  diffuses  away  as  it  passes  down¬ 
stream.  Further,  what  is  required  for  a  description  of  electromagnetic 
scattering  (and  to  a  lesser  extend  for  the  prediction  of  self-luminous  radia¬ 
tion)  is  the  description  in  local  space  of  the  distributions  of  the  time  de¬ 
pendent  properties  in  the  frequency  domain  and  their  behavior  with  free  stream 
conditions. 

Tatar  ski,  V.I. ,  Wave  Propagation  in  a  Turbulent  Medium,  McGraw-Hill,  1961. 
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z.  THE  PHENOMENOLOGY  OF  WAKE  TRANSITION 


The  question  of  stability  of  flows  to  small  pe rturbations  is  not  the 
question  of  transition.  However,  the  inherent  instability  of  laminar  mo¬ 
tion  at  sufficiently  high  Reynolds  numbers,  R,  is  responsible  for  the  ultimate 
transition  to  turbulent  flow.  The  value  of  the  transition  Reynolds  num¬ 
ber,  depends  not  only  on  the  minimum  critical  Reynolds  number,  Rcr  . 
determined  from  stability  theory  for  infinitesimal  disturbances*  but  also 
on  the  initial  size  of  the  actual  disturbances  with  the  most  '’dangerous” 
frequencies. 

The  approach  adopted  here  to  the  problem  of  the  inherent  insta¬ 
bility  of  steady  fluid  dynamic  flows  and  the  onset  of  turbulence  is  the 
classical  one.  Conceptually,  if  small  perturbations  which  inevitably 
occur  in  the  flow  tend  to  increase  with  time,  then  the  flow  is  absolutely 

unstable.  The  absolute  instability  of  the  flow  for  R  >  R  leads  to  the 

c  r 

appearance  of  a  nonsteady  flow  with  characteristic  period  and  with  small 
but  finite  amplitude  which  increases  with  R,  This  periodic  flow  usually 
can  be  described  in  terms  of  a  succession  of  vortex  filaments.  When  the 
Reynolds  number  increases  further,  a  time  comes  when  the  periodic  flow 
in  turn  becomes  unstable  with  the  appearance  of  a  new  period.  When  the 
Reynolds  number  increases  still  further,  more  and  more  new  periods 
appear  in  succession.  The  range  of  Reynolds  numbers  between  succes¬ 
sive  appearances  of  new  frequencies  diminishes  rapidly  in  size;  the  new 
flows  themselves  are  on  a  smaller  and  smaller  scale.  The  flow  rapidly 
becomes  complicated  and  confused  and  we  say  that  such  a  flow  is  turb¬ 
ulent. 

Consider  steady  wake  flow  in  the  x -direction.  A  slight  perturba¬ 
tion  whose  wave  components  in  the  course  of  time  are  amplified  is 
carried  downstream.  The  energy  in  the  velocity  fluctuations  builds  up 
in  the  amplification  regime  into  the  ordered  time  dependent  vortex  flow 
with  characteristic  frequencies.  Over  a  small  range  of  Reynolds  numbers 
this  vortex  system  may  persist  to  infinity,  decaying  through  molecular 
diffusion  alone.  With  subsequent  increase  in  R,  two  things  happen:  (1) 

The  effective  origin  —  the  downstream  position  at  which  there  is  a 
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certain  level  of  energy  density  at  the  ordered  frequency  (and  its  har- 

>  s' 

iconics)  —  oi  the  vortex  system  moves  upstream,  towards  the  body. 

{ Z )  the  velocity'  fluctuation  wake  energy  in  the  ordered  flow  ie  continu- 
-aiiy  passed  down,  with  increasing  distance  downstream .  to  the  energy 
in  the  random  small  scale  eddy  velocity  fluctuations  containing  all  fre¬ 
quencies,  'Finally',  an.  B.  •"*  01 ,  .  the  boundary  layer  on.  the  body  is  fully  __ 
turbulent  and  the  distribution  of  fluctuation  energy  density  in  the  base  region 
already  has  finite  values  at  all  frequencies. 

In  a  comprehensive  cet  of  experiments  at  low  speeds  (rearward 

3  4  5 

facing  step-,  rectangular  jet,  and  flat-plate  wake,  )  Sato  has  shown  that 

transition  in  a  free  shear  layer  essentially  follows  the  hierarchy  de¬ 
scribed  above:  initiation  by  Tollmein  waves  that  grow  exponentially  in 
the  downstream  direction,  followed  by  a  nonlinear  regime  characterised 
by  a  double  row  of  vortices,  and  finally  a  region  in  which  the  fluctuation 
loses  regularity  and  gradually  develops  into  turbulence. 

Leonard  and  Keck's  modification  of  the  conventional  schlieren 
system  uses  a  sodium  vapor  in  the  test  medium  with  an  interference 
filter  in  the  optical  domain  with  peak  transmission  20  R  off  the  5896  & 
sodium  line.  ^  The  extremely  large  value  of  the  index  of  refraction  at 
the  resonance  line  greatly  enhances  the  effective  index  of  the  test  medium 
for  the  light  passing  through  the  filter.  At  the  specified  conditions-,  it  is 
at  least  two  orders  of  magnitude  larger  than  that  for  air  in  the  optical 
domain.  This  makes  it  possible  to  extend  the  pressure  regime  for  good 
sfcclieren  contrast  from  that  of  a  few  percent  of  atmospheric  to  pressures 
less  than  half  of  one  percent,  and  to  considerably  increase  the  detail  apparent 
in  pictures  taken  in  the  two  to  three  percent  atmospheric  region.  Fur¬ 
ther,  using  a  constant  partial  pressure  of  sodium  vapor,  the  sensitivity 
of  the  schlieren  system  is  constant  and  thus  independent  of  the  gas 
pressure. 

Sodium  schlieren  pictures^  allow  the  observation  of  the  develop¬ 
ment  of  structure  in  the  hypersonic  sphere  wake  at  the  onset  of  transi¬ 
tion.  It  is  hypothesized  that  qualitatively  the  development  is  similar  to 
the  classical  one  described  above.  Near  the  body,  small  perturbations 
in  the  flow  are  in  the  form  of  Tollmien-Sehlichting-like  wa”es  which,  in 
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the  course  of  time,  are  amplified  and  carried  downstream.  These  waves 
are  longitudinal  in.  orientation  and  of  several  body  diameters  in  length. 
The  fluctuating  wake  energy  builds  up  in  the  succeeding  nonlinear  ampli¬ 
fication  regime,  and  this  growth  process  leads  into  a  nearly  ordered 
flow  at  a  set  of  characteristic  frequencies.  This  structure  might  Taede- 
scribed  as  a  succession  of  vortex  filaments,  which  in  the  sphere  case 
is  like  a  procession  of  vortex  loops,  one  tied  to  another. Proceed¬ 
ing  downstream,  there  is  a  continual  degradation  of  die  fluctuating  wake 
energy  from  the  nearly  ordered  nonsteady  flow,  with  significant  energy 
at  distinct  frequencies,  down  into  random  small  scale  fluctuations  con¬ 
taining  all  frequencies  and  proceeding  towards  fully  developed  turbulence. 
Figure  1  is  a  composite  of  five  separate  sodium  schlieren  pic¬ 
tures,  each  single  picture  taken  on  a  separate  run.  All  runs  were  made 
at  constant  fixed  range  pressure  of  1.7  cm  Hg  of  Fox  the  body 

size  used  in  these  experiments,  d  -  0.  56  cm,  this  pressure  is  in  the 

a 

range  of  the  onset  of  hypersonic  wake  transition  to  noneteady  flow. 
Approximately  0,375  mm  Hg  partial  pressure  of  sodium  vapor  was  main¬ 
tained  in  the  test  section.  The  pellet  velocities  varied  expe rimen  tally  by 
about  plus  or  minus  10  percent  from  the  mean  of  14,000  ft/sec.  (Hyper- 

l 

sonic  wake  transition  has  been  shown  to  be  insensitive  to  10  percent 

8  9 

changes  in  velocity.  *  7)  The  spark  source  was  triggered  by  a  variable 
time  delay  after  the  passage  of  the  projectile,  thus  providing  observa- 

i 

tions  of  various  portions  of  the  wake  as  shown.  The  knife  edge  was 
horizontal,  thus  emphasizing  vertical  gradients. 

The  photographs  in  Fig.  1  illustrate  the  hierarchy  of  transition 
modes.  Horizontal  wave-like  striations  of  several  body  lengths  long  are 
plainly  visible  originating  near  to  or  from  the  base  region.  These  cor¬ 
respond  to  the  small  perturbations  mentioned  above.  The  second  picture 
shows  these  perturbations  after  some  fluid  mechanical  amplifications. 

The  striations  are  now  formed  into  a  nearly  ordered  flow  and  the  struc¬ 
ture  appears  to  oe  like  a  procession  oi  vortex  filaments  having  loop-like 
shape.  Weak  shocks  are  visible  emanating  from  some  of  these  filaments, 
indicating  their  large  scale  and  their  large  velocity  relative  to  the  in- 
viscid  flow.  In  the  third  photograph,  the  original  perturbations  have  been 


mm 


Sodium  schlieren  composite  of  a  0.22-inch  sphere  hypersonic 
wake  in  nitrogen.  Pressure  =  1.7  cm  Hg,  V  ~  14,000  ft/se< 


amplified  to  a  large  degree  and  the  pattern  looks  more  or  less  like 
vortex  loops,  one  stacked  behind  the  other,  “exElbitifig  t&eliearly  periodic 
component  to  toe  flow.  The  next  picture  in  the  sequence  shows  some 
indication  of  the  deterioration  of  the  large  scale  vortex  filaments  as 
they  fill  up  with  email  scale  eddies.  ~The  IajBtZ.pcture^ 

nation  of  this  decay  process  towards  fyliy  developed  turbulence.  The . 

shock  waves  seen  in  the  last  three  pictures  are  Tefl.ee ti'ons  of  the  pro- 
je..  tile's  bow  shock  from  the  test  section  walls.  This  set  of  pictures 
□ear  transition  illustrates  the  origin,  growth,  and  the  beginning  of  the 
breakdown  towards  fully  developed  turbulence  of  hypersonic  nonsteady 
wake  flow.  _  . 

Other  sensitive  schlieren  systems  omitting  the  sodium  vapor  mhd- 

ification  have  also  produced  photographs  of  hypersonic  wakes  and  each 

data  will  be  used  in  the  following  sections.  At  the  Reynolds  number 

corresponding  to  the  data  shown  in  Fig.  1,  such  schlieren  photographs 

appear  to  show  a  .  rapid  development  of  turbulence  downstream  of  the 

a 

body  at  a  point  30  body  diemeters  from  the  pellet.  Figure  1  suggests 

*  i  -''l*"-’-" . 

that  what  is  seen  at  30  body  diamete r a  dowmsfcream  (second  photograph) 
is  the  continued  simplification  oi  the  small  density  gradient  perturbations 
already  present  near  the  body  in  the  zero  to  10  body  diameter  region, 
as  shown  in  the  first  photograph.  It  is  suggested  that  the  increased 
refractivity  of  the  sodium  resonance  radiation  has  allowed  us  to  detect 
weak  density  fluctuations  where  they  first  appear  right  near  to  the  body. 
The  threshold  of  the  air  schlieren  systems  used  in  other  experiments 
may  not  have  allowed  observation  of  the  original  small  fluctuations  until 
after  some  fluid  mechanical  amplification.  Only  when  the  density  gradi¬ 
ents  become  sufficiently  large  to  create  a  light  fluctuation  in  such  an 
air  system  will  turbulent  structure  appear  in  such  a  photograph.  Advanced 

schlieren  techniques  now  confirm  the  gradual  development  of  turbulence 
41 

downstream. 

Figure  2  shows  a  photograph  of  the  self-luminous  wake  (taken  by  the 
racetrack  technique*®)  at  a  similar  shoulder  Reynolds  number  to  the  sodium 
schlieren  results  of  Fig.  1.  The  luminosity  of  the  near  wake  region  appears 
to  have  a  long  wavy  motion  which  downstream  at  about  30  body  diameters 


develops  into  large  scale,  vortex  filaments.  This  amplification  of  the  near 
wake  wave  motion  into  large  scale,  vortex  filaments  shown  by  the  racetrack 
technique  in  Fig.  2,  is  presumed  to  be  the  same  process  as  the  amplification 
of  the  small  perturbations  near  to  the  body  into  vortex  filaments  evident  by 
the  sodium  schlierea  technique  it,  “Fig:  i'«  •Th  this-sense.  -bo&r  observations" 
support  the  thesis  that  this  amplification  does  not  portray^ a  sudden  transition 
to  turbulent  flow.  ’  8La^ei,~thfix.e-ia’a-  hierarchy  of  Jfiow  patterns  beginning  an 
the  body  flow  region  and  continuing  downstream  which  lead  to  the  fully  devel¬ 
oped  turbulent  wake  structure. 

Figure  3  shows  a  high  quality  shadowgraph  behind  a  cone.  The 
shadowgraph  technique  is  sensitive  to  the  second  derivative  of  density 
and  hence  empha  sizes  smaller  scale.  Also,  in  the  cone  flow  field,  the 
entropy  changes  from  the  boundary  layer  -  induced  inner  wake  to  the 
shock-induced  outer  wake  are  greater  than  in  the  sphere  case.:  -Thus , 
shadowgraphs  of  cone  wakes  in  the  transition  regime  may  be  expected 
to  show  more  detail  of  the  development  of  fluid  mechanical  unsteadiness 
than  schlieren  photographs  of  sphere  wakes.  ^  One  sees  in  the  details 
of  wake  transition  in  Fig,  3  the  same  sort  hierarchy  of  unsteady  flow 
regimes  as  described  above  in  experiments  with  other  diagnostics. 

The  hypersonic  wake  has  three  strong  compressibility  effects 
which  must  be  considered  in  addition  to  the  shape  of  the  velocity  profilet 
pressure  gradient,  Mach  number,  and  lateral  temperature  variation. 
ToUmein*^  has  shown  that  the  existence  of  a  point  of  inflection  in  the 
velocity  profile  is  a  necessary  and  sufficient  condition  for  the  amplification 
of  infitisimal  disturbances.  In  boundary  layer  flow  with  favorable  pressure 
gradients,  the  velocity  profile  has  negative  curvature  everywhere.  Thus, 
these  boundary  layer  flows  have  higher  transition  Reynolds  numbers  than 
those  with  adverse  pressure  gradients  since  the  latter  possess  velocity 
profiles  with  inflection  points. 

13 

Wake  velocity  profiles  have  inflection  points  a  priori.  Pretsch 
has  shown  that  the  influence  of  the  pressure  gradient  enters  stability 
analysis  only  through  the  velocity  profile,  u(y).  Thus,  in  considering  the 
stability  of  wake  flow  with  pressure  gradient,  one  need  only  consider  the 
pressure  effect  on  the  distribution  of  the  velocity  profile  with  y. 
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For  slender  cones,  like  total  cone  angle  les3  than  3 
show  that  the  p  -*■  point  is  close  to  the  end.  of  tl 


_  u  -  c 

<1  (3) 

e 

for  the  development  of  instabilities,  where  c  is  the  phase  velocity  of 

s 

the  disturbance  in  body  coordinates.  Numerical  calculations  by  Gold 
indicate  that  ior  hypersonic  flow  this  is  satisfied  not  too  far  downstream 
of  the' body  (<  lOOd)  in  both  the  blunt  and  slender  body  cases. 

18 

It  should  be  pointed  out  that  in  recent  work  by  Mack  on  the  two- 
dimensional  problem,  even  when  no  neutral  or  self-excited  subsonic 
disturbances  corresponding  to  the  fundamental  mode  exist,  he  found  the 
existence  of  self-excited  disturbances  of  the  higher  modes.  Energy  trans¬ 
fer  for  these  modes  may  be  provided  by  the  '’bouncing1'  of  the  waves 
back  and  forth  within  the  supersonic  region. 

14 

The  effect  of  a  hot  wake  core  on  the  stability  analysis  is  two-fold: 
(1)  As  the  wake  core  temperature  excess  increases,  the  range  of  rela¬ 
tive  Mach  numbers  over  which  unstable  disturbances  can  exist  also 

increases;  (2)  however,  if  the  relative  Mach  number  is  less  than  M  , 

,  cr 

the  amplification  rates  per  wavelength  for  two-dimensional  flows  of  a 
hot  wake  decrease  significantly  with  increasing  AT  and  are  rather  unaf¬ 
fected  in  the  axisymmetric  case. 

For  a  blunt  body,  the  condition  in  Eq.  (3)  is  satisfied  close  to  the 

re  compression  region  well  before  the  p  -*  p^  position  is  reached.  Thus, 

the  end  of  the  compressibility  effects  in  hypersonic  wake  transition  for 

blunt  bodies  is  controlled  by  the  pressure  decay.  For  a  slender  cone 

one  set  of  calculations  show  that  at  high  altitudes,  say  from  100,000  to 

200,000  ft,  >1  may  persist  for  some  several  hundred  body  diam- 

8  s  19 

eters  downstream.  Thus,  the  end  of  the  compressibility  effects  in 
hypersonic  wake  transition  for  cones  is  controlled  by  the  M  decay. 

Since  the  compressibility  effects  are  related  to  the  hypersonic 
processes  nea.r  to  the  body,  the  part  of  the  wake  which  is  dominated  by 
these' local  details  will  be  called  the  "Near  Wake.  "  By  definition  the 
"Far  Wake"  will  be  where  the  wake  velocity  profiles  have  reached  their 
asymptotic  behavior,  property  gradients  have  become  small  and  details 
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of  .body  shape  have  been  lost.  The  Interpolation  regime  will  be  the  regime 
where  one  traverses  from  the  Near  to  the  Far  wake:  where  the  strong 
favorable  pressure  gradient,  the  strong  density  gradients,  and  the  effect 
of  body  shape  all  become  less  important  as  one  proceeds  downstream. 

Before  proceeding  with  analysis  of  the  wake  transition  experiments, 
we  discuss  those  cases  where  there  has  been  injection  of  mass  from  the 
body  into  the  wake.  Such  injection  could  be  done  deliberately,  such  as  would 
occur  in  some  types  of  boundary  layer  control  methods,  or  could  be  the  re¬ 
sult  of  surface  ablation.  No  laboratory  experiments  on  this  effect  are  avail¬ 
able  in  the  literature. 

We  can  make  certain  remarks  on  what  might  be  expected  in  the  case 
of  low  mass  injection  rates  into  the  flow  field.  We  will  assume,  based  on 
certain  peripheral  experimental  information  to  be  presented,  that  if  the  rate 
of  mass  injection  is  small  relative  to  the  rate  of  mass  of  air  encountered  by 
the  body,  that  there  is  little  ablation  effect  on  wake  transtion.  This  assump¬ 
tion  has  some  theoretical  base.  Consider  an  experiment  in  wind  tunnel  co¬ 
ordinates.  It  is  reasonable  to  assume  that  the  injected  mass  leaves  the  body 
surface  at  small  velocity  relative  to  the  fast  free  stream.  Thus,  the  momen¬ 
tum  balance  for  the  control  volume  surrounding  the  body  requires  that  all  of 
the  momentum  gained  by  the  injected  mass  in  the  free  stream  direction 
be  lost  by  the  free  stream,  which  maintains  the  total  drag  of  the  body 
constant.  (We  neglect  for  this  argument  the  effect  of  ablation  on  the 
boundary  layer  profile  which  will  affect  the  skin  friction  drag.  )  Although 
the  momentum  defect  passing  out  of  the  control  volume  (which  is  equal 
to  the  body  drag)  is  unchanged,  the  velocity  distribution  may  be  locally 
altered  so  as  to  produce  a  more  unstable  profile.  However,  as  Gold** 
has  calculated,  small  differences  in  profile  shape  have  small  effects  on 
growth  of  disturbances.  Thus,  the  assumption  of  small  effects  when  the 
ratio  of  injected  mass  rate  to  intercepted  is  small  seems  to  be  reason¬ 
able.  Further  remarks  will  be  made  at  appropriate  places  in  the  following 
discussion  for  the  effects  to  be  expected  in  the  case  of  high  mass  injection 
rates. 
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3..  THE  spectral  outline  for  the  experimental  data 

The  spectral  intensity  of  the  wake  fluctuations  may  be  roughly 
divided  between  the  discrete  part  at  the  characteristic  frequencies 
(ordered)  and  the  random  part  at  all  other  frequencies  (disordered). 

After  the  minimum  critical  Reynolds  number  has  been  passed,  increas¬ 
ing  the  Reynolds  number  moves  upstream  the  position  at  which  some 
given  distribution  of  wake  fluctuations  between  ordered  and  disordered 
is  found.  Figure  4  illustrates  this  process  schematically.  For  the 
present  purpose  of  illustrative  discussion,  pressure  decay,  Mach  num¬ 
ber  and  temperature  effects  are  neglected.  The  vertical  axis  represents 
the  distribution  of  powt:  >  between  ordered  and  random  (disordered)  mo¬ 
tion.  The  two  axes  in  the  horizontal  plane  are:  (1)  a  representative 
Reynolds  number,  R^;  and  (2,'  the  distance  downstream,  x/d.  The  up¬ 
stream  progression  is  seen  to  be  the  locus  in  the  horizontal  plane  cor¬ 
responding  to  a  given  power  spectral  distribution  (ordered-disordered). 
Lines  obtained  by  planar  cuts  for  the  two  other  processes  are  also  shown 
the  degradation  of  wake  fluctuations  down  from  the  nearly  ordered  flow 
proceeding  towards  fully  developed  turbulence:  (1)  at  fixed  x/d  with  in¬ 
creasing  R^;  and  (2)  at  fixed  R^  with  increasing  x/d.  If  pressure  decay, 
Mach  number,  temperature  and  ablation  effects  are  important,  there  is 
a  separate  surface  in  (R^,  x/d,  spectral)  space  for  each  combination  of 
values  for  these  four  effects. 

It  is  useful  to  think  of  the  experimental  techniques  employed  for 
hypersonic  wake  transition  work  as  mapping  out  different  horizontal 
contours  of  the  (R^,  x/d,  ordered-disordered)  surface.  Close  to  the 
value  (1.0  -  0.0)  on  the  ordered-disordered  spectral  scale  are  tech¬ 
niques  which  are  sensitive  to  property  fluctuations  or  powers  thereof. 

For  example,  self-luminous  radiation  is  sensitive  to  a  high  power  of 
temperature.  Very  slight  fluctuations  in  temperature  will  produce  per¬ 
ceptible  structure  in  a  self-luminous  photograph.  The  appearance  of 

structure  in  the  self-luminous  photographs  with  increasing  R,  seems  to 

8  ^ 

have  the  characteristics  of  an  "off-on"  event.  fech  would  be  the  nature 
of  the  phenomenon  represented  by  the  horizontal  locus  corresponding  to 
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Schematic  illustration  of  development  of  turbulence 


a  spectral  valu>  just  below  (1.  0-0.  0).  Further  down  the  spectral  scale 
would  be  the  in  e rfe romete r  technique  which  is  directly  related  to  density 
fluctuations. 

Next  down  the  spectral  scale  is  the  schlieren  technique.  The 
schlieren  technique  is  sensitive  to  the  first  derivative  of  the  density. 
Thus,  the  schlieren  system  registers  unsteadiness  when  the  density 
gradients  (which  are  associated  with  the  smaller  scale  eddies)  become 
sufficiently  sharp  to  create  light  fluctuations  which  are  above  the  thres¬ 
hold  level  to  produce  structure  in  the  schlieren  photograph. 

At  still  lower  values  on  the  spectral  scale,  comes  the  shadow¬ 
graph.  The  shadowgraph  is  sensitive  to  the  second  derivative  of  density 
and  hence  emphasizes  still  smaller  scale  eddies  than  the  schlieren 
technique. 

If  these  fcur  optical  techniques  were  used  on  the  same  experiment, 
each  would  map  out  the  loci  of  points  corresponding  to  four  horizontal 
contours  in  Fig.  4.  Unfortunately,  it  is  not  presently  possible  to  obtain 
quantitative  information  about  spectral  distributions  from  all  of  these  tech¬ 
niques.  However,  their  relative  locations  with  respect  to  each  other  on 
the  spectral  axis  should  be  in  the  order  in  which  they  were  discussed 
above. 

Diagnostics  in  other  parts  of  the  electromagnetic  spectrum  are 
also  used;  e.g.,  microwaves,  specially  for  hypersonic  wake  transition. 
Microwave  data  are  fox  electron  wake  unsteadiness.  It  is  most  likely 
that  fluid  mechanical  unsteadiness  is  a  prerequisite  for  electron  density 
unsteadiness.  Thus,  it  is  to  be  expected  that  fluid  mechanical  unsteadi¬ 
ness  is  as  close  or  closer  to  the  body  than  microwave  data  indicates. 
Further,  the  selected  location  of  the  beginning  of  fluid  mechanical  un¬ 
steadiness  determined  from  visual  examination  of  photographs  is  not 
sensitive  to  the  wavelength  of  the  unsteadiness  at  that  point.  However, 
the  location  of  wake  unsteadiness  determined  by  microwave  scatter¬ 
ing  is  very  sensitive  to  the  scale  of  the  electron  density  unsteadiness 
as  represented  by  the  refractive  index  fluctuations.  The  wake  after 
having  gone  fluid  mechanically  unsteady  may  not  produce  the  smaller 
electron  fluctuation  scale  for  microwave  scattering  until  further  down- 
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stream  or  until  a  higher  ambient  pressure  is  reached.  Thus,  at  a  given 
pressure,  the  relative  location  of  microwave  wake  transition  further 
downstream  than  the  fluid  mechanical  observation  is  a  reasonable  state 
of  affairs,  even  without  considering  the  problem  of  electron  production. 
The  results  of  ballistic  range  experiments  at  General  Motors  Defense 
Research  Laboratory^®  support  this  argument.  Figure  5  shows  this  rela¬ 
tionship  for  both  sphere  and  cone  experiments;  that  is,  the  microwave 
transition  curve  is  further  downstream  or  at  higher  ambient  pressure 

than  the  curve  traced  out  by  the  schlieren  techniques. 

21 

Demetriades  used  the  integrated  rms  output  of  a  hot-wire  as  a 
diagnostic  in  hypersonic  cylinder  wake  experiments.  For  x^,  he  arbi¬ 
trarily  selected  the  location  of  either  the  enhancement  or  the  maximum 
in  the  rms  output.  The  (R^.x/d)  locus  traced  by  the  criterion  of  the 
enhancement  of  the  rms  signal  would  correspond  to  a  higher  position  on 
the  spectral  axis  than  the  criterion  of  the  maximum  output.  Hypersonic 
cylinder  wake  transition  is  different  from  that  for  three-dimensional 
bodies.  In  the  cylinder  case,  time  to  pass  information  along  the  wake 
shear  layer  in  the  direction  perpendicular  to  the  two-dimensional  planar 
flow  is  long  compared  with  time  to  pass  information  across  the  layer 
and  with  time  to  flow  downstream  around  the  body.  This  is  not  the  case 
for  wake  shear  layers  of  bodies  of  revolution.  Thus,  the  details  of  cylin¬ 
der  hypersonic  wake  transition  are  different  from  those  in  the  axisym- 
metric  case  and  the  cylinder  case  is  not  treated  here.  The  concepts 
used  in  this  paper  to  analyze  hypersonic  wake  transition  for  bodies  of 
revolution  were  actually  developed  first  for  cylinder  wake  transition.  The 
results  are  qualitatively,  though  not  quantitatively,  the  same.* 

To  locate  the  (Rd#x/d)  locus  at  its  actual  spectral  value,  one 

would  have  to  resolve  the  fluctuation  spectrum.  Roshko  did  this  with 

22 

hot  wire  techniques  in  low  speed  cylinder  wake  experiments.  So  far 
this  has  not  been  done  in  hypersonic  experiments,  although  a  start  has 

"5  Goldburg ,  A.  ,  Unpublished  prepared  comment  presented  on 

Demetriades'  1963  paper21,  AIAA  Hypersonics  Conference,  Cambridge, 
Mass.  ,  August  1963. 
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Fig.  5  Comparison  of  schlieren  and  radar  wake  transition  location 
for  both  spheres  and  cones.  (General  Motors  Defense 
Research  Laboratory,  Ref.  20,  Fig.  10) 
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been  made  by  Kendall  at  JPL23  and  Hromas  at  TRW.  *  However,  it  is  pos¬ 
sible  to  qualitatively  relate  the  results  of  Roshko  to  the  data  obtained  from 
schlieren  diagnostics  in  ballistic  ranges. 

Slattery  and  Clay9  carefully  analyzed  their  schlieren  photographs  of 
hypersonic  wakes  in  the  transition  regime.  For  each  ballastic  range  shot 
they  measured  the  distance,  xg,  from  the  body  to  the  point  where,  to  their 
schlieren-eye-mind  analysis  system,  some  perceptible  amount  of  structure 
could  be  seen  in  the  photograph;  e.  g.  ,  see  Fig.  3.  Figure  6  includes  a  pres¬ 
entation  of  their  data  for  spheres  at  fts  7.  5.  The  ordinate  is  xg/d  and 
the  abscissa  is  shoulder  ’  the  ^eyn°l<is  number  calculated  for  conditions 
corresponding  to  those  on  the  centerline  streamline  after  it  goes  through  the 
shock,  around  the  nose  and  is  turned  back  to  horizontal.  For  blunt  bodies, 
this  is  90°  from  the  stagnation  point  vertical  and  for  cones  this  is  one-half 
the  total  cone  angle  from  the  end  of  the  cone  surface.  The  shoulder  flow  con¬ 
ditions  are  taken  as  typical  of  the  near  wake  flow  conditions.  The  local  Rey¬ 
nolds  number  per  unit  length  along  the  edge  of  the  viscous  core  in  the 
near  wake  is  more  or  less  constant  and  is  near  the  shoulder  value  for 
blunt  bodies  and  equal  to  that  for  slender  cones.  Since  the  shoulder 
values  are  easier  to  calculate  and  correspond  to  a  definite  position  in 
space,  the  shoulder  conditions  are  chosen  for  convenience  and  for  cer¬ 
tain  reasons  which  will  become  apparent  later  on. 

xg  is  considered  in  this  paper  in  the  sense  of  being  that  down¬ 
stream  distance  required  for  the  development  of  unsteadiness  into  that 
degree  of  turbulence  with  sufficient  sharpness  of  density  gradient  fluc¬ 
tuations  to  register  above  the  threshold  level  of  the  diagnostic  technique. 

It  is  assumed  that  this  degree  of  turbulence  corresponds  for  each  sepa¬ 
rate  reading  (each  separate  run)  to  some  unknown  but  more  or  less 
constant  value  on  the  spectral  axis.  Fig.  4.  Thus,  the  data  in  Fig.  6 
is  looked  upon  as  mapping  out  some  particular  horizontal  contour  of  the 
(Rd,xg/d,  spectral)  surface:  f(Rd,xg/d)  =  0. 

From  the  quantitative  hot-wire  fluctuation  spectral  data  given  by 
22 

Roshko,  it  is  possible  to  draw  the  curves  for  his  low  speed  experi¬ 
ment  for  the  80  percent  -  20  percent  (orde red-disordered)  plane  and 
the  20  percent-80  percent  (ordered-disordered)  plane.  It  is  seen  in  Fig.  6 

*  Rrivate  communication,  L.  Hromas,  September,  1966. 
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Comparison  of  x8/d,  data  for  hypersonic^  and  low-speed 
experiments. 22  Curves  for  Roshko's  incompressible  cylinder 
results  are  from  his  hot-wire  spectral  data.  8 


that  the  shape  of  the  hypersonic  transition  curve  (Xg/d.R^)  is  exactly 
similar  to  the  loci  of  points  traced  out  by  the  results  of  the  quantitative 
hot-wire  data  in  the  low  speed  experiments  of  Roshko.  Thus,  we  are 
encouraged  that  the  spectral  framework  for  the  data  as  illustrated  sche¬ 
matically  in  Fig.  4  is  appropriate  for  an  analysis  of  the  hypersonic  wake 
transition  experiments. 

4.  THE  CONCEPTUAL  OUTLINE  FOR  THE  SCALING  LAWS 

Hypersonic  wake  transition  derives  its  separateness  from  low 
speed  wake  transition  from  the  fact  that  in  hypersonic  wake  transition 
there  are  two  wakes  left  by  the  body:  the  shock-induced  wake  and  the 
boundary  layer-induced  wake.  Fig.  7.  The  author  prefers  this  nomen¬ 
clature  due  to  its  specific  nature,  clearly  labelling  the  origin  of  each 
vortical  flow  and  the  mechanism  for  drag  deposition  in  the  flow  field. 
Other  names  in  use  respectively  are:  outer-inner,  entropy -viscous ;  and 
inviscid-viscous.  The  author  feels  that  the  last  two  are  clearly  misnomers 
and  are  definitely  misleading  to  a  large  number  of  interested  technologists, 
in  that  both  wakes  are  "entropy"  and  "viscous"  wakes,  and  that  the  shock- 
induced  wake  is  clearly  not  without  viscosity.  The  "entropy"  wake  ‘  is 
so-called  because  for  blunt  bodies  most  of  the  shock  drag  is  produced 
by  the  part  of  the  shock  within  the  sonic  line  which  is  where  the  large 
entropy  change  takes  place  across  the  shock.  However,  in  the  case  of 
slender  bodies  most  of  the  entropy  change  takes  place  in  the  boundary 
layer  and  so  this  nomenclature  is  not  appropriate.  The  "inviscid"  wake 
is  so-called  because  in  calculating  normal  forces  on  bodies  in  supersonic 
flow,  it  is  customary  to  separate  the  flow  field  into  a  "free  stream",  in 
which  viscous  effects  are  dropped  as  unnecessary  for  the  normal  surface 
force  problem,  and  into  a. boundary  layer  where  viscous  effects  are  used  to 
satisfy  the  tangential  no- slip  condition  at  the  surface.  However,  the  wake 
transition  problem  is  a  free  shear  layer  problem  in  a  real  gas  possessing 
finite  viscosity  involving  a  transition  Reynolds  number,  and  should  not  in 
any  sense  be  spoken  of  or  thought  of  as  being  "inviscid",  without  viscosity. 

There  is  a  part  of  small  disturbance  theory  associated  with  the  Orr- 
Sommerfeld  equation  in  which  the  viscous  terms  are  dropped.  24  The 


RECOMPRESSION  SHOCK 


SHOCK  INDUCED  WAKE 

BNtfY  LAYER  induced  wake 


Flow  field  nomenclature.  Shock -induced  wake  (SiW)  and 
boundary  laye  r -induced  wake  (BLiW). 


7 


mathematical  problem  of  an  "inviscid”  stability  calculation  must  be  rec¬ 
ognized  as  associated  with  laminar  flow  which  are  solutions  of  the  full 
Navier -Stokes  equations  where  the  influence  of  viscosity  on  the  mean 
flow  has  been  taken  into  account.  What  is  neglected  in  the  inviscid 
infinitesimal  small  disturbance  theory  is  the  influence  of  friction  on 
the  infinitesimal  fluctuation.  Of  course,  as  soon  as  the  fluctuation 
builds  up  into  finite  size  in  the  non-linear  amplification  regime  going 
towards  turbulent  transition,  viscosity  plays  a  role.  Therefore,  in 
discussing  transition,  the  use  of  the  '‘inviscid’*  nomenclature  is  unfortu¬ 
nate  and  will  not  be  used  here. 

With  two  contiguous  wakes,  the  shock  induced  wake  (SiW)  and  the 
boundary  layer  induced  wake  (BLiW),  there  are  three  transition  possi¬ 
bilities:  1)  the  SiW  has  gone  unsteady  and  the  BLiW  has  not;  2)  the  BLiW 
has  gone  unsteady  and  the  SaW  has  not;  3)  both  wakes  have  gone  unsteady 
either  in  unison  or  in  separate  motion.  It  is  possible  that  when  one 
wake  has  gone  unsteady  it  may  "trip”  the  other.  For  example.  Fig.  8, 
a  self-luminous  racetrack  photograph,  shows  both  wakes  evidently  in  un¬ 
steady  sinuous  motion,  each  independent  of  the  other.  By  using  xenon  as 
the  test  gas,  very  high  free  stream  Mach  numbers  (M^  «  25)  were  ob¬ 
tained,  and  hence  very  high  flow  field  temperatures  were  produced  behind 
the  bow  shock.  This  produces  gas  luminosity  in  the  SiW  as  well  as  the 
stronger  luminosity  produced  by  the  vaporized  ablation  products  in  the 
BLiW,  thus  visually  distinguishing  between  the  two  wakes.  Evidence  in 
support  of  the  possibility  of  the  SiW  being  in  unsteady  motion  independ¬ 
ently  of  the  BLiW  was  first  given  by  Fay  and  Goldburg  (Ref.  8, 

Section  E). 

With  the  above  nomenclature  and  with  the  delineation  of  the  wake 
by  definition  into  three  parts:  Near  Wake,  Interpolation,  and  Far 
Wake,  we  next  build  a  conceptual  outline  for  the  experimental  data. 

The  Reynolds  number  will  be  adopted  as  the  correlating  non- 

dimensional  parameter  for  the  wake  transition  studies:  This  is  taken 

directly  from  all  previous  fluid  mechanical  studies  and  experience.  It  is 

the  primary  oarameter  derived  from  the  , non-dimen sionalization  of  the 

24 

Navier-Stokes  equations  of  motion,  and  has  had  great  usefulness  in  the 
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Fig.  8  Luminous  undulating  hypersonic  wake  with  the  SiW  and  the 
BLiW  in  apparently  uncorrelated  time -dependent  motions. 
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study  of  transition  in  boundary  layers,  jets  and  incompressible  wakes. 
The  numerical  value  of  the  Reynolds  number  is  determined  by  the  local 
properties  characteristic  of  the  flow  in  the  vicinity  of  transition. 

In  the  Far  Wake,  the  fluid  has  relaxed  its  memory  of  the  details 
of  the  way  in  which  momentum  was  lost  to  the  body  (wind  tunnel  body 
fixed  coordinates  will  be  used  unless  otherwise  stated).  The  velocity, 
density  and  temperature  gradients  have  become  small;  the  property  dis¬ 
tributions  have  reached  their  asymptotic  shape.  For  laminar  flow,  the 
BLiW  by  molecular  diffusion  has  spread  out  and  smeared  itself  into  the 
SiW  so  that  the  profiles  show  no  sharp  gradients  within  the  interior  of 


the  total  wake. 

In  the  Far  Wake,  it  will  be  postulated  that  the  whole  wake  goes 
unsteady  together  at  some  unknown  station  or  series  of  stations  way, 
way  downstream.  This  will  be  the  first  onset  of  wake  unsteadiness, 

~  00  ,  and  the  total  shear  layer  representing  the  whole  momentum  defect 
equal  to  the  total  drag  of  the  body  will  participate  in  the  transition.  Thus, 
the  appropriate  length  scale  will  be  the  total  ^ C^A  of  the  body,  or  equiv¬ 
alently  the  total  momentum  thickness  of  the  total  wake,  9,  being  a  sum 
of  that  of  the  SiW  and  the  BLiW.  Thus,  the  Reynolds  number  for  tran¬ 
sition  to  unsteadiness  in  the  Far  Wake  will  be  based  on  ^  CpA: 


r  x  =■  constant.  Figure  9a  presents  this  discussion  schemati- 

callyv  The  shaded  area,  represents  the  diffusion  of  the  BLiW  into  the 
SiW,  thus  effectively  making  one  wake  with  total  momentum  defect  equal 
to  the  total  drag  of  the  body. 

In  the  Hear  Wake,  the  two  wakes  are  clearly  delineated  by 
sharp  property  gradients  across  the  edge  of  the  BLiW.  The  two 
wakes  have  strong  memory  of  iheir  past  origins  and  are  strongly  influ¬ 
enced  by  the  local  details  of  the  flow  field  around  the  body.  It  has  been 
observed  that  to  schlieren  diagnostics  the  SiW  usually  appears  smooth 
in  this  region,  while  the  BLiW  has  gone  unsteady.  In  the  schematic  this 
state  of  affairs  is  represented  in  Fig.  9b.  It  can  be  seen  in  such  a 
case  that  it  may  be  possible  for  the  BLiW  to  be  unsteady  and  for  the 
SiW  at  very  high  entropies  to  be  opaque  to  the  diagnostic.  In  such  a 
case  as  this,  the  diagnostic  would  not  be  able  to  follow  the  upstream 
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movement  of  the  x  beyond  the  point  where  the  unsteadiness  enters 
within  the  high  entropy  region  of  the  SiW,  which  remains  steady.  In 
other  words,  if  the  steady  entropy  region  of  the  SiW  is  opaque  to  the 
diagnostic,  the  unsteady  BLiW  would  not  be  apparent  until  it  broke 
through  the  "entropy"  line.  Laboratory  evidence  for  this  phenomenon 
was  first  obtained  by  Wilson.  ^ 

In  the  Near  Wake,  local  hypersonic  phenomena  depending  on  body 

shape  dominate  the  flow  field  as  is  the  case  for  all  body  related  shear 

27 

layers  —  boundary  layers  and  separation  shear  layers.  Hidalgo,  and 

21 

Professor  Lees  and  his  graduate  students  extended  the  transition  cor¬ 
relation  ideas  developed  for  body  related  shear  layers  and  proposed 
that  the  wake  transition  correlation  length  would  be  the  distance  from 
the  origin  of  the  Near  Wake  to  the  transition  point.  For  convenience, 
the  origin  for  x^  is  taken  from  the  base  of  the  body.  The  end  of  the 
recompressicn  region  might  be  more  appropriate,  but  since  values  of 
x  of  interest  are  usually  greater  than  two  body  diameters  downstream, 
and  since  the  correlations  to  be  developed  will  be  valid  to  no  better 
than  a  factor  of  two,  the  choice  of  origin  is  not  critical  since  the  re- 

28 

compression  region  itself  is  like  one  or  two  body  diameters  downstream. 

Further,  oy  continuing  analogy  they  suggested  that  R  based  on  xtf  would 

be  a  strong  function  of  the  local  Mach  number  at  the  edge  of  the  BLiW 

at  the  location  of  x^.M^.  Thus,  equals  a  constant,  Cj^,  which 

varies  with  M  .  tr 

e 

In  the  Interpolation  regime,  one  traverses  from  the  Near  to  the  Far 
Wake.  Going  downstream,  the  strong  favorable  pressure  gradient,  the  strong 
density  gradients,  the  large  Mach  number  and  temperature  changes  across  the 
boundary  layer  and  shock -induced  wakes  decrease  and  become  less  important. 
Going  upstream  in  this  region,  the  wakes  begin  to  take  on  their  separate  char¬ 
acters  and  the  unsteadiness  becomes  more  and  more  to  be  located  mostly  in 
the  BLiW.  This  situation  is  represented  schematically  in  Fig.  9c. 

In  the  Interpolation  regime,  the  Navier-Stokes  equations  will  be  used 
with  certain  simple  arguments  to  derive  an  interpolation  curve  to  go 

from  the  R  — r-  correlation  of  the  Far  Wake  to  the  Rv  correlation 

J  CDA  xtr 

of  the  Near  Wake.  The  simplest  such  interpolation  formula  is 
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(c“A« j  *  I1  •  IM/^  I] 

(  )  j,  is  the  transition  value  for  the  onset  of  Far  Wake  transition  and 

a  is  some  fitting  parameter-  to  be  dele r  mined  by  the  data,___ Thus,  con¬ 
sider  a  given  body  in  a  given  flow  field*.  For  otxset  of  wake  unsteadiness 
far  downstream  in  the  Far  Wake  (Fig.  9a)  x.  —*  *®  ,  the  Far  Wake  law 
is  obtained 

53  t  R^Jc . a  )  *  constant,  {b) 

For  R  per  unit  length  increasing  to  very  large  values,  x^  decreasing, 
the  transition  point  moving  towards  the  body,  the  Near  Wake  law  is 
obtained. 


At  a  given  flight  velocity  for  gas  mixtures  gove-ned  by  binary 
collisions,  it  is  possible  to  represent  the  Reynolds  number  as  the  free 
stream  pressure  times  the  characteristic  length.  Thus,  the  hypersonic 
wake  transition  laws  become 


Far  Wake:  Vs?  =  constant 


lNear  Wake:  =  constant  (function  of  M  (8) 

M  ,  the  Local  Mach  number  at  the  edge  of  the  BLdW  is  here  replaced 
c 

by  tne  Mach  number  at  the  sho ’older  of  the  body  after  the  flow  has  been 
turned  back  to  the  horizontal,  see  sketch  in  Fig.  10.  Since  in  the 

near  wake  region,  does  not  vary  sigrificantly  as  a  function  of  dis¬ 
tance  downstream  along  the  BLiW  edge,  and  since  the  local  Mach  number 
at  the  shoulder  of  the  body  i3  close  :o  M  ,  M  ,  is  selected  here  as  be- 
ing  sufficiently  characteristic  of  the  local  hypersonic  flow  field  in  the 
region  near  the  particular  shaped  body. 

The  correlations  (7)  and  (8)  lend  themselves  to  presentation  on  a 
log -log  plot  with  x^r  as  ordinate  and  as  abscissae.  Fig.  10.  Thus, 
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Section  7. 


5.  THE  FAR  WAKE  CORRELATION 
The  Far  Wake  is  that  region  where  all  details  of  the  body  shape 
have  diffused  away.  Only  the  SiW  is  discernible,  the  BLiW  having 
diffused  into  it.  Fig.  9;  all  that  remains  is  information  about  the  total 
momentum  defect  deposited  by  the  body  ~  C^A.  Transition  in  the  Far 
Wake  is  here  characterized  as  an  off-on  phenomenon.  As  indicated  in 
Fig.  10  we  wish  to  find  for  each  free  stream  pressure,  pM  ,  a  size  of 
body  whose  wake  will  show  unsteadiness  somewhere  very  far  downstream 
as  that  pressure  threshold  is  passed.  In  this  step,  no  attention  is  paid 
to  the  hierarchy  of  the  transition  modes.  A  scaling  law  is  sought  for 
the  onset  of  wake  time  dependency  in  terms  only  of  the  total  momentum 
defect  deposited  by  the  body  drag. 

29  l - 

Results  from  incompressible  flow  indicate  that  the  totai-IC^A 

left  by  the  body  is  the  correlating  length  parameter  independent  of  the 

means  of  deposition;  i.  e.  ,  regardless  of  the  details  of  the  combination 

of  normal  pressure,  skin  friction,  base  and  induced  pressure  drags, 

Arkhipov^  performed  the  small  disturbance  stability  analysis  for  a 

1  ? 

planar  wake  with  velocity  profiles  taken  from  the  Tollmein’s  solutions.’’ 


The  results  of  his  wake  stability  analysis  show  that  the  parameter  which 
characterizes  a  transition  of  the  flow  to  a  state  of  instability  is  Reynolds 
number  based  on  total  wake  momentum  thickness:  RQ  53  24,  where  mo¬ 
mentum  thickness  is  related  to  C_A  by 


P«ui  JrM+1  =  Drag  =  l/Zp*  uj  CDA 


j  -  Q  =  planar 
j  =  1  =  circular 


The  suggestion  that  the  onset  of  wake  unsteadiness  is  characterized 
by  a  transition  value  of  R^»  independent  of  the  mechanism  of  drag  pro¬ 
duction,  is  also  supported  by  experimental  work  on  cylinders  and  flat 
plates.  For  cylinders,  wake  transition  occurs  at  Rd  *  50,  with  corre¬ 
sponding  Rg  «  42  (Ref.  25).  At  this  value  of  Red,  the  total  drag  is 
divided  approximately  half  and  half  between  pressure  (form)  drag  and 
skin  friction  drag  (Ref.  25,  Fig.  154).  For  flat  plates  of  length  l , 
Taneda's  experiments^*  give  wake  transition  at  Rj  8  1000,  with  corre¬ 
sponding  R  <=*44.  In  the  case  of  the  flat  plate  the  total  drag  is  produced 

v 

by  skin  friction  alone.  The  fact  that  Arkhipov's  theoretical  value  is  lower 
than  the  experimental  values  is  to  be  expected  since  the  experiments  deal 
with  transition  to  finite  unsteadiness  whereas  the  theory  deals  with  the 
onset  of  amplification  of  infinitesimal  disturbances.  Another  considers- 
ation  is  that  the  theoretical  work  assumes  for  the  velocity  profile  an 
asymptotic  Gaussian  form;  whereas  in  the  experiment,  the  instabilities 
may  develop  nearer  to  the  flat  plate  where  the  velocity  profile  is  not 
Gaussian.  The  actual  velocity  profile  at  the  onset  of  amplification  may 
give  a  higher  theoretical  Rfi  transition  than  the  result  developed  from  the 
assumed  Gaussian  shape. 

29 

In  Goldburg  and  Florsheim,  transition  in  incompressible  wakes 
of  three-dimensional  bodies  was  considered:  spheres  and  slender  cones 
of  various  cone  angles.  They  showed  for  the  three-dimensional  case  that 
although  the  (Rd)  T  for  the  bodies  tested  varied  by  a  factor  of  two,  when 
the  total  momentum  thickness  or  was  used  as  the  length  param¬ 

eter  for  onset  of  wake  unsteadiness:  (R^  =  95  was  a  correlation  to 
better  than  10  Dercent. 


Thus,  there  is  theoretical  and  experimental  support  for  correlating 
the  onset  of  hypersonic  Far  Wake  unsteadiness  with  the  Reynolds  number 
based  on  a  length  parameter  characteristic  of  the  total  drag  area  of  the 
body  represented  by  the  total  momentum  defect  across  the  whole  SiW, 
including  the  diffused  BLiW.  {This  correlation  was  first  used  for  hyper¬ 
sonic  wakes  by  W.  Schultis  in  an  npublished  report:  Study  of  transition 
from  laminar  to  turbulent  flow  in  wakes,  Nov.  1961.  )  The  ballistic 
range  laboratory  data  for  all  ballistic  range  experiments  available  in  the 
open  literature  is  plotted  in  Fig.  11  to  determine  empirically  the  con¬ 
stant  in  Eq.  (7),  pM  =  constant. 

Consider  first  the  drum  camera  photographs  of  self-luminous  wakesj 
These  data  characte ristically  show  either  a  smooth  diffuse  texture  or 
show  gradients  in  the  luminosity  (Fig.  12).  As  z.  function  of  increasing 
free  stream  pressure,  luminous  structure  is  either  absent  or  then  pres¬ 
ent.  No  upstream  motion  of  structure  with  increasing  Reynolds  number 
is  usually  seen.  The  absence  or  presence  of  luminous  structure  with 
increasing  Reynolds  number  seems  to  have  the  character  of  an  off-on 
event,  as  if  the  transition  "snapped-in"  from  very  far  downstream  almost 
instantaneously.  Such  would  be  the  character  of  the  phenomenon  repre¬ 
sented  by  the  horizontal  locus  corresponding  to  a  spectral  value  just 
below  {1.0  -  0.0)  on  the  orde red-disordered  spectral  scale  in  Fig.  4. 

It  is  suggested  that  self-luminous  radiation  structure  in  the  near 
wake  is  a  sensitive  test  for  onset  of  finite  transition  in  the  Far  Wake. 

The  reason  for  this  suggestion  is  that  radiation  goes  as  a  very  high 
power  of  the  temperature.  Any  small  anomalies  in  temperature  will  give 
rise  to  greatly  amplified  structure  in  the  radiation  pattern.  Thus,  very 
small  flow  field  anomalies  in  the  Near  Wake,  which  will  be  amplified 
into  the  onset  of  finite  unsteadiness  in  the  Far  Wake,  may  be  detectable 
in  the  Near  Wake  with  self-luminous  radiation.  Figure  10  sets  forth  the 
transition  information  in  air  for  the  self-luminous  diagnostic  on  a 
p-  and  plot.  Laminar  corresponds  to  a  self-luminous 

photograph  without  structure;  turbulent,  with  structure;  and  intermediate 
means  that  the  observer  could  not  make  a  definite  determination  because 
of  some  very  slight  indications  of  unsteadiness. 
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Fig.  12 


R  cetrack  photographs  of  luminous  hypersonic  wakes, 

d  --  0. 22".  . 
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The  schlieren  values  for  transition  were  obtained  by  the  investi¬ 
gator  laying  a  marker  on  the  photograph  where  perceptible  unsteady 
structure  in  the  wake  was  apparent.  Fig.  3.  These  data  ail  showed  the 
behavior  as  represented  in  Fig*  6*  Thus,  as  inauiejit  pressure  was 
decreased  from  run  io  run,  the  or  xg  marked  by  the  investigator 
moved  further  and  further  back  in  th,  wake  until  some  minimum  pres¬ 
sure  was  passed  below  which  no  structure  was  ever  apparent  in  the  wake. 
The  schlieren  data  are  represented  in  Fig.  11  as  a  vertical  line  at  the 
approximate  velocity  of  the  set  of  experiments.  At  the  top  of  the  line  are 
the  initials  of  the  performing  laboratory.  At  the  bottom  is  S  or  C,  for 
sphere  or  cone  data.  Beside  the  line  is  the  number  of  runs  represented. 
The  top  of  the  line  is  terminated  with  an  arrow  at  the  highest  value  of 
P«  Vc  qA  for  that  set  of  experiments.  The  arrow  indicates  that  the  wake 
would  continue  turbulent  with  increasing  values  of  Obviously, 

this  limit  could  be  pushed  higher  and  higher,  eventually  producing  tran¬ 
sition  in  the  body  boundary  layer.  If  a  lower  asymptote  is  reached  for 
a  minimum  pressure  (Reynolds  No. )  for  transition,  such  as  shown  in 
Fig.  6,  the  arrow  is  terminated  at  the  minimum  p^  Vc£a  of  the  ex¬ 
periment  set  with  a  horizontal  mark.  For  example,  the  line  at  8000  ft/ 
sec  represents  the  data  in  Fig.  6.  If  no  minimum  asymptote  appears 
on  such  a  plot  as  Fig.  6,  the  line  is  terminated  at  the  bottom  with  a 
downward  facing  arrow  indicating  that  the  experiment  set  did  not  pro¬ 
ceed  to  low  enough  values  of  pw  to  reach  the  disappearance  of 

wake  unsteadiness. 

Figure  11  includes  all  available  laboratory  ballistic  range  data. 
Transition  is  a  statistical  phenomenon  and  the  off-on  characterization 
used  for  the  Far  Wake  lends  itself  to  a  probability  analysis.  The  tran¬ 
sition  results  of  Fig.  11  give 


MEAN 

1  a 

2  a 


-  1.0  cm  x  cm  Hg 
=  +  0.  5 

=  +  0.75 


We  choose  here  to  take  the  low  la  value  as  the  constant  for  the  fhr 

Wake  t  ran  sit  ion  correlation  of  F*  ,  |7), 


-35  - 


P  co 


0.  5  cm  Hg  x  cm 


(9) 


It  is  emphasised  that  the  yjC^A  in  the  Far  Wake  correlation  (9)  is 

the  total  yjC^  A  and  corresponds  to  all  of  the  drag  in  the  SiW  and  in 

the  BLiW,  the  latter  having  diffused  in  some  fluid  mechanical  sense 

into  the  form  :r  in  the  Far  Wake  such  that  they  both  go  unsteady  as  a 

whole  or,  if  still  discernible,  together.  This  correlation  is  meant  to 

include  angle  of  attack  effects  as  well,  so  long  as  total  yJC^A  at  the 

angle  of  attack  is  used  in  the  correlation. 

The  use  of  ablating  bodies  {the  self-luminous  photographs  depend 

on  the  luminosity  of  the  ablation  products  for  the  observable  diagnostic) 

for  transition  work  is  questionable.  However,  in  these  experiments  with 

small  ablation  rates  (Calculations  show  that  for  these  experiments  the 

ratio  of  the  mass  rate  of  ablation  with  respect  to  the  mass  rate  of  air 

encountered  is  less  than  0(10  ^).  )  the  transition  value  of  p^ 

does  not  vary  between  individual  sets  of  experiments  with  and  without 

ablating  bodies.  If  ablation  rates  were  large,  like  of  the  order  of 

the  mass  rate  of  air  encountered,  one  might  expect  (1)  either  an 

alteration  to  he  characteristic  length  used  (perhaps  including  some 

2 

characterization  of  the  "force",  m  ,  .  U  ,  divided  by  l/2  p  U  to 

^  abl  c0  r  °°  «> 

make  an  addi  oonal  "length"  )  or  (2)  an  alteration  which  is  a  function 

of  the  ablation  rate  to  the  value  of  the  constant  in  the  Far  Wake  law. 

26 

Qualitative  experiments  by  Wilson  support  the  above  arguments. 

Since  free  shear  layers  are  extremely  unstable  and  have  critical 
Reynolds  numbers  of  the  order  of  10,  rarefaction  effects  may  be  im¬ 
portant.  Representing  the  Reynolds  number  by  Mach  number  times  the 
inverse  of  tin  Knudsen  number,  we  see  that  for  high  Mach  number  flows, 
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say  like  thi  order  of  10,  that  the  Knudsen  number  is  like  order  one. 
Thus,  hype  sonic  wake  unsteadiness  occurs  very  close  to  the  end  of  the 
continuum  flow  regime,  if  not  somewhat  into  the  rarefied  flow  regime. 
The  scaling  constant  given  in  Eq.  (9),  based  on  ballistic  range  data, 
assumed  for  spheres  and  cones  the  continuum  value  of  Cj-j.  It  may  be 

that  C~>  C„  .  at  the  onset  of  wake  unsteadiness  by  10  or 

D  D  continuum  7 

20  percent. 

In  conclusion,  it  is  pointed  out  that  the  correlation  constant  of 

Eq.  (9)  is  rood  only  to  like  a  factor  of  two.  This  is  readily  seen  in 

Fig.  11  where  the  transition  region  is  shown  shaded  to  a  factor  of  two 

either  side  of  the  mean.  For  hypersonic  transition  work  this  is  quite 

usual.  Fc.  example,  the  transition  Reynolds  number  for  heat  transfer 

to  blunt  no  ed  bodies  in  hypersonic  flow  is  also  known  only  to  a  factor 
33 

of  two.  ■',/hat  one  is  dealing  with  in  compressible  turbulent  transition 
is  a  random,  spectrally  disorganized  highly  nonlinear  process,  the 
onset  of  which  is  decidedly  influenced  by  the  anomalies  in  the  steady- 
state  condition  of  the  ambient  environment. 

6.  THE  NEAR  WAKE  CORRELATION 

In  the.  Near  Wake,  local  hypersonic  phenomena  depending  on  body 
shape  dominate  the  flow  field.  The  strong  favorable  pressure  gradient, 
the  large  temperature  and  Mach  number  differential  across  the  SiW 
tend  to  stabilize  it.  Usually  only  the  BLiW  is  observed  to  go  into  finite 
unsteadiness  in  schlieren  and  shadowgraph  photographs.  The  SiW  is 
like  a  smooth  sheath  around  the  corrugated  BLiW  and  for  the  electro¬ 
magnetic  scattering  observable  could  specularly  reflect  the  signal  and 
thus  shield  the  BLiW  from  the  radar  diagnostic. 

As  explained  in  Section  3,  the  point  at  which  the  investigator  lays 
down  the  marker  on  the  schlieren  photograph  to  mark  the  onset  of 
finite  unsteadiness  in  the  BLiW  is  considered  here  in  the  sense  of  being 
that  downstream  distance  required  for  the  development  of  disorder  in 
the  flow  with  that  degree  of  turbulence  which  possesses  sufficient 
sharpness  of  density  gradient  fluctuations  to  register  abewj  the  three- 
hold  level  of  the  diagnostic  technique.  It  is  assumed  that  this  degree 
of  turbulence  corresponds  for  each  separate  reading  (each  separate  run 
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of  each  sep.  rate  investigator)  to  some  unknown  but  more  or  less  con¬ 
stant  value  on  the  spectral  axis.  Fig.  4.  Thus,  we  assume  that  the 
Near  Wake  .-chlieren  data  are  information  on  how  some  level  of  non-discrete 
fluctuating  power  moves  upstream  with  increasing  Reynolds  number. 

It  is  common  for  transition  in  Mach  number  stabilized  body-related 
free  shear  ayers  that  the  hierarchy  of  fluid  mechanical  states  from 
laminar  to  ally  developed  turbulence  is  compressed  into  a  very  small 
streamwise  distance,  of  the  order  of  two  or  three  shear  layer  thick¬ 
nesses.  This  corresponds  to  the  sharp  drop  off  in  percentage  ordered 
motion  with  x/d  at  high  in  Fig.  4.  This  is  in  contradistinction  to 
Far  Wake  transition,  where  the  streamwise  distance  for  transition  may 
be  of  the  order  of  several  hundred  shear  layer  thicknesses.  This  corre¬ 
sponds  to  the  gentle  drop-off  in  ordered  motion  with  x/d  at  low  in 
Fig.  4.  The  fact  that  most  flow  visualization  data  of  the  wake  transi¬ 
tion  process  is  in  the  near  wake  has  misled  many  people  into  believing 
that  transition  is  an  a  priori  precipitous  process.  The  whole  concept  of 
for  wakes  is  largely  due  to  an  overemphasis  on  near  wake  transi¬ 
tion  picture  u. 

27  21 

Hidalgo  at  Avco  Everett  first  suggested  and  Demetriades  and  Gold 

first  demonstrated  conclusively  that  the  Rxtr  ~  constant  correlation  employed 

in  transition  of  body  related  shear  layers  could  be  carried  over  to  the  near 

wake.  The  s.rong  dependence  of  this  correlation  on  local  wake  Mach  number, 

27 

M  ,  at  the  point  of  wake  transition  {also  suggested  by  Hidalgo  ),  was  dem- 

C  or  |Q  o  j  o  / 

or.strated  by  Demetriades  and  Gold,  Zeiberg,  and  Pallone.  Others  ’  * 

have  indicated  a  dependence  on  Mn  for  cones,  but  for  cones  this  is  essentially 

M  .  A  correlation  on  M^  does  not  characterize  the  essence  of  the  problem 

in  the  local  wake  flow  field.  Hence,  it  does  not  correlate  both  spheres  and 

cones  where  the  local  Me  is  quite  different  for  the  two  geometries;  in  the 

case  of  spheres,  M  ,  not  being  close  to  Ma, 

In  the  Near  Wake  region,  calculations  show  that  Mg  does  not  vary 

significantly  as  a  function  of  distance  downstream  along  the  wake  edge, 

and  th  f  the  local  Mach  number  at  the  shoulder  of  the  body  (see 

sketch  ir  ’0  for  definition  of  )  *8  c*ose  tr  -*e*  Thero,.<re. 

here,  M  ,  .  .  is  selected  as  a  parameter  related  to  body  shape  wh. _n 

shoulder  r 


-38- 


is  sufficient  y  characteristic  of  the  local  hypersonic  flow  field  in  the 
region  near  the  particular  body.  M  .is  high  and  close  to  M  for 
slender  cones,  and  is  low  and  close  to  3  for  spheres. 

The  x  data,  as  determined  and  published  by  the  investigators 

1  9  19  3a 

from  Lincoln  Laboratory,  Avco/RAD,  Naval  Ordnance  Laboratory, 

and  General  Motors^  are  now  used  on  a  p  •  x  vs  M  ,  plot,  p  •  x 

r°°  tr  sh  r  tr 

represents  through  binary  scaling  the  Ry  .  The  correlation  shown  in 
Fig.  13  is  *r 


p  •  x. 
r°°  tr 


,  J1  +  .  082  M  ,  )  .  „  . 

=  10  sh  [  cm  Ilg  •  cm  ] 


This  determines  the  constant,  a  function  of  Mach  number,  for  the  Near 

Wake  Law  of  Eq.  (8).  It  is  noted  that  the  Near  Wake  correlation  of  (10) 

is  also  good  to  like  a  factor  of  two.  Thus,  to  each  value  of  M  ,  .  ,  , 

°  shoulder 

one  can  assign  a  line  of  negative  one  slope  on  the  (x  ,  )  plot,  Fig.  10. 

It  should  be  pointed  out  that  in  the  Near  Wake  correlation,  care 

was  exercised  not  to  use  data  points  in  the  so-called  "sticking"  re- 

26 ' 38  ■  ■ 

gion.  ’  There  is  an  upstream  limit  of  the  Near  Wake  transition 
correlation  where  the  recompression  region  interacts  with  the  wake 
transition  pr  ocesses  to  retard  further  upstream  motion  of  x^^  independ¬ 
ent  of  furthe  r  increases  in  pressure.  This  clearly  is  not  included  in  an 

Rv  correlation, 
tr 

Now,  conceptually,  the(xtr,  )  plot.  Fig.  10,  can  be  filled  with 

a  series  of  vertical  lines  assigning  a  (pw  for  each  and  a 

series  of  negative  one-slope  lines  assigning  an  (pM  .  x^)  for  each 

M  ,  ,  ,  .  What  is  needed  finally  is  an  interpolation  curve  to  connect 
shoulder  _ 

the  far  wake  onset  altitude  for  the  yJC^A  of  the  body  with  the  near 
wake  Rxtr  behavior  for  the  of  the  body. 

7.  THE  INTERPOLATION  REGIME 

In  the  Interpolation  Regime,  one  traverses  from  the  Near  to  Far 
Wake.  The  strong  favorable  pressure  gradient,  the  strong  density  gra¬ 
dients  and  the  effect  of  body  snape  all  become  less  important.  Thus, 
compressibility  effects  have  smaller  and  smaller  intxuences  in  this 
regime,  and  we  look  at  the  unsteady  incompressible  N' vie  r -Stokes 
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equation  for  guidance.  Time  is  nondimensioaalized  with  the  time  re¬ 
quired  for  ths  decaf  of  the  nearly-ordered  wake  flow  down  into  some 
given  spectral  distribution  of  random  small-scale  turbulence,  tb  =  xc/Uoo 
and  is  the  distance  to  a  given  spectral  distribution  of  unsteadiness, 
here  corresponding  to  The  size  of  the  large  scale  eddy,  which  is 

of  the  order  of  the  size  of  the  shear  layer,  is  taken  as  the  character¬ 
istic  length  in  the  y -direction.  These  big  eddies  do  not  grow  appre¬ 
ciably  by  molecular  diffusion  in  the  time  required  to  decay  through  growth 

g 

of  instabilities.  The  size  of  the  large  scale  is  like  d  for  blunt  bodies 
i  39 

and  like  d/3  lor  slender  cones.  The  distance  required  for  the  fluctua¬ 
tion  level  to  build  up  to  a  specified  finite  amplitude,  say  A  a,  is  taken 
as  the  characteristic  length  in  the  x -direction,  and  here  is  equivalent 
to  x  . 

The  linearized  time  dependent  equation  of  motion  in  the  x-direction 
is: 


■> 

d  "  u' 


9  u‘  .  -  d  u*  ,  .  3u  _  ,a“u  , 

~  ft  +  *  -y  -  *  <77  — 


UD 


,vhe  re 


u  =  u 


with  the  usual  notation  for  mean  velocity,  u,  and  time  dependent  vel¬ 
ocity.  u*.  v'  is  of  the  order  of  u',  u1  =  a(x)  •  f(t)  and  r  is  the  time 
required  for  *h"  fluctuation  level,  to  build  up  to  an  amplitude  Aa. 
Substituting  into  Eq.  (11): 
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whe  re 


denotes  the  nondimensionalized  functions  and  their  deriva- 


I  1 


tives. 


Wake  correlation  for  increasing:  ~  a  constant  dependent 

on  Mach  number.  Thus,  for  the  purpose  of  the  hypersonic  wake,  Eq.  (15) 
is  modified  into 


(Jb±  \ 

A  -  1R,M 

[  k  i 

V  *tr/ 

i  JX  «  g 

\  &  / 

or  with  binary  scaling  and  with  d  replaced  by 


(  CM  \ 

u  /,  .  (p-  V 

fCSX  >T  \ 

\  p«  xtr  ) 

1  \  P.  yJCZX  ) 

This  relationship  is  now  fully  determinate  (the  a  of  Eq.  (4)  having 
come  out  to  be  unity)  since  we  have  determined  the  correlation  constants: 

as  given  on  the  right  hand  side  of  correlation  (10)  and  the  value  of 
<P»  )  as  given  on  the  right  hand  side  of  correlation  (9).  Consid¬ 

ering  the  data  for  a  single  body,  Eq.  (18)  can  be  simplified  to 

x,  . 

tr  1 


P OQ  'PoO 


The  dimensional  analysis  above  having  been  presented,  we  take  the  par¬ 
ameters  of  Eq.  (19)  as  the  variables  of  a  relationship  to  be  determined 
empirically  by  the  transition  data  itself: 


tr 

CM,/p°°T 


e  =  P«/PcoT 


or  or.  a  log-log  plot 


r|'  =  log  r|;  p’  =  log  p 

It  is  seen  that  the  two  asymptotes  are  correctly  given: 

Far  Wake:  V  =  00  ,  p 1  =  0:  p  -  (p^)^ 

for  the  correct  of  the  experimental  vehicle. 


Near  Wake: 


P'  =  0: 


x  =  C, 
tr  1 


the  correct  C.,  corresponding  to  M  ,  of  the  experimental  vehicle,  Thus, 
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the  experimental  data  correlated  on  the  variables  t)1  and  p*  will  de¬ 
termine  a  single  normalized  interpolation  curve  for  all  intersections  of 
the  Near  and  Par  Wake  correlations.  1.  e«  ,  for  every  vehicle  in  every  flow 
field  each  realization  of  the  x^,  curve  is  characterized  by  one  y  C^A 
and  one  —  Fig.  10. 

The  x  data  in  the  Interpolation  regime  as  determined  and  puh- 
t*  9  /  19 

listed  by  the  investigators  from  Lincoln  Laboratory,  Avco/RAD, 

_  _  _  _ *■»  j _ _ _ *■»  /■  _ _ _ 

Naval  Ordnance  Laboratory,  and  General  Motors6  are  now  used  on 
a  rj',  p*  plot.  Fig.  14.  The  reason  there  is  so  little  cone  data  for  the 
Interpolation  regime  as  well  as  for  the  Far  Wake,  Fig.  11,  is  that  gun 
launched  cones  fly  straight  and  level  at  high  pressures.  There  is  an 
obvious  sparsity  of  data  for  transition  onset  in  the  downstream  part  of 
the  wake,  which  requires  experiments  at  lower  pressure.  Most  transi¬ 
tion  data  have  been  collected  for  the  Near  Wake  regime  because  cones 
tend  to  be  dynamically  stable  at  high  ambient  pressures  and  schiieren 
definition  is  good.  Almost  all  of  the  published  cone  data  are  for  the 
Near  Wake  transition  onset,  Fig.  13.  In  Fig.  14  these  data  line  up 
along  the  Near  Wake  asymptote:  p'  +  rj'  =  0,  and  are  not  shown. 

The  Interpolation  regime  data  in  Fig.  14  are  used  to  establish  an 
Interpolation  curve  to  connect  the  Far  Wake  law,  (7)  and  (9),  and  the 
Near  Wake  law,  (8)  and  (10);  see  Fig.  10.  The  Interpolation  curve 
shown  in  Fig.  14  is  normalized  for  all  intersections  of  the  Far  Wake 
correlation  (on  the  size  of  vehicle,  yC^A  )  and  the  Near  Wake  corre¬ 
lation  (on  the  shape  vehicle  in  any  flow  field,  M^).  Thus,  it  is  now 
possible  to  construct  a  singie  smooth  curve  which  includes  the  repre¬ 
sentation  of  the  Far  Wake  correlation  and  the  representation  of  the  Near 
Wake  correlation  and  which  has  a  shape  not  too  dissimilar  from  Eq.  (19) 
and  the  form  shown  in  Fig.  6  for  Roshko's  data. 

Of  course,  we  do  not  know  what  is  the  quantitative  nature  of  the 
spectral  distribution  whose  forward  movement  is  being  represented  by 
these  correlations.  We  have  assumed  that  each  time  an  investigator 
lays  down  a  marker  on  the  schiieren  photograph  for  the  onset  of  wake 
unsteadiness,  that  each  event  corresponds  to  approximately  the  same 
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spectral  stimuli.  Lake  the  other  correlations,  the  Interpolation  curve 
of  Fig.  14  is  good  to  like  a  factor  of  two. 

Recently,  investigations  at  Lincoln  Laboratory41  and  General 
Motors  '  have  drawn  attention  to  phenomena  which  calls  into  question 
the  simplicity  and  the  smoothness  of  the  proposed  transition  curve  start¬ 
ing  ^11  the  way  from  the  Far  Waite  to  and  through  the  Interpolation  re¬ 
gime  into  and  terminating  at  the  end  of  the  Near  Wake  region;  as  rep¬ 
resented  in  Fig.  14.  To  assess  the  nature  of  these  considerations,  we 
plot  all  of  the  sphere  data  in  Fig.  15.  We  use  the  given  axes  in  order 
to  include  data  which  were  presented  by  Slattery  and  Clay,  only  as  a 
set  of  curvesonthis  particular  grid  (1962,  Fig.  5).  Clay,  Labitt  and 
Slattery  claim  that  their  high  Mach  number  data  (M^  =  18-20)  pro¬ 
duces  a  separate  curve  from  some  of  their  data  at  lower  Mach  number 
(Ref.  41,  Fig.  6).  However,  if  all  of  their  lower  Mach  number  data 
and  if  all  of  the  General  Motors  data  is  included  in  Fig.  15,  there 
docs  not  show  two  or  several  separate  curves  but  just  one  smooth, 
broad  bandof  data,  with  spread  like  a  factor  of  two,  with  shape  repre¬ 
sented  by  the  Near-Interpolation-Far  Wake  curve  as  developed  in  the 
present  paper.  In  fact,  it  was  identically  these  data  which  were  used 

to  establish  the  correlations  developed  herein. 

41 

Clay,  Labitt  and  Slattery  call  attention  to  a  second  feature  of 
their  data.  When  they  consider  their  high  speed  data  alone,  there  ap¬ 
pears  to  be  a  sharp  break  to  the  left  at  about  2000  ball  diameters. 

26 

Wilson  believes  his  data  also  evidence  such  behavior  (1965,  Fig.  7). 
However,  we  see  in  Fig.  15  that  when  we  place  ail  of  the  sphere  data 
or.  one  grid,  ail  that  remains  is  a  smooth  band  of  data  following  the 
behavior  of  the  correlations  and  interpolation  developed  here.  As  more 
and  more  data  are  collected,  one  would  expect  this  band  to  become  solid 
with  data  pc’ots. 

A  third  effect  was  noticed  by  Wilson  (1965)26  for  very  high  Mach 
number  tests,  Mm  >  20.  For  blunt  bodies  at  high  Mach  numbers, 

there  is  a  region  of  high  entropy  production  across  the  nearly  normal 
part  of  the  shock,  Fig.  7.  This  region  called  the  entropy  wake  is  one 
of  high  thermal  excitation  and  resulting  low  density.  The  sensitivity  of 
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Fig,  15  All  ballistic  Range  Sphere  Transition  Data 


schiieren  systems  is  proportional  to  the  density;  and  as  a  consequence  , 
when  the  entropy  wake  ts  hot  enough,  as  it  is  at  very  high  >  20, 

then  it  becomes  “opaque"  to  the  schiieren  diagnostic  and  shields  the 
turbulence  in  the  BLiW  from  schiieren  observation.  Under  these  ex¬ 
perimental  conditions  —  blunt  bodies  at  high  Mach  number  with  a  schiieren 
diagnostic  —  the  turbulent  BLiW  would  not  become  apparent  until  it  had 
"broken  through"  the  low  density  entropy  regime  of  the  SiW.  Thus,  even 


though  the  x  had  moved  upstream  in  the  BLiW,  the  shielding  of  the 

I  i 


schiieren  diagnostic  by  the  entropy  wake  would  make  the  apparently 


remain  fixed  at  the  point  where  BLiW  entered  the  high  entropy  region. 
Wilson’s  "break-through"  data  in  Fig.  15  illustrates  this  phenomenon. 
However,  it  is  important  to  realize  that  this  phenomenon  is  clearly  un¬ 
derstandable  in  terms  of  the  shielding  of  the  BLiW  by  the  SiW  which 
hides  from  the  observer  the  continuing  forward  motion  of  fluid  mechan¬ 
ical  transition. 


S.  THE  HYPERSONIC  WAKE  TRANSITION  MAP 


We  now  construct  quantitatively  the  hypersonic  wake  transition  map 
which  was  shown  schematically  in  Fig.  10.  The  transition  map  will  pro¬ 
vide  a  method  for  making  the  transition  curve,  Ktr»  as  a  function  of 
pressure  ratio  (altitude),  for  a  given  body  size  —  V^A,  and  a  given 
body  shape  at  a  given  flight  velocity  — 

The  Far  Wake  transition  law  was  determined  to  be 


=  2.  2  x  IQ'4  ft 


(9) 


Then,  to  each  altitv.de,  there  wiil  correspond  a  body  size,  , 

from  Eq.  (9),  such  that  the  'sake  of  that  body  when  travelling  at  hyper¬ 
sonic  speeds  down  through  that  attitude  will  go  unsteady.  On  the  Altitude- 
xtr  Map  of  Fig.  10,  this  is  represented  as  a  vertical  line  at  each  alti¬ 
tude  coming  down  from  xtr - -  00  for  the  appropriate  ylC^A.  This  has 

been  d^r.e  in  Fig.  lb. 

The  Near  Wake  transition  law  was  determined  to  be 
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Fig.  16  Hypersonic  Wake  Transition  Map. 
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where  x  is  in  centimeters.  Hence,  on  the  (p*,  ixtr)  Map  of  Fig.  16 
there  are  a  series  of  negative  one-slope  lines,  each  for  a  different 
To  facilitate  the  use  of  the  Transition  Map,  there  is  given  a 
conversion  chart  and  graph  for  to  Mg^,  Fig.  17. 

,  The  Interpolation  curve  developed  in  the  previous  section.  Fig.  14, 
is  also  shown  in  the  map,  Fig.  16.  For  ease  of  use,  the  Interpolation 
curve  is  included  with  this  report  as  a  clear  plastic  overlay. 

9.  CONCLUSIONS 

The  Transition  Map,  Fig.  16,  provides  a  method  for  making  the 
transition  curve,  x^  as  function  of  pressure  ratio  (altitude),  for  any 
given  body  at  any  given  flight  velocity:  Connect  the  Far  Wake  onset 
altitude  for  ihe  yjC^A  of  the  body  with  the  Near  Wake,  p^  •  xtr,  behavior 
for  the  corresponding  Mgh  of  the  body.  Thus,  Fig.  16  becomes  some¬ 
thing  of  a  complete  Wake  Transition  Map  for  a  large  range  of  Mach 
numbers  and  body  sizes  and  body  shapes. 

This  map  is  based  on  an  analysis  which  takes  the  Reynolds  num¬ 
ber  aa  the  correlating  parameter  and  divides  the  wake  transition  into 
three  regions:  1)  Far  Wake,  where  details  of  body  shape  have  diffused 
away  and  the  characteristic  length  is  the  total  of  the  vehicle; 

2)  Near  Wake,  where  local  hypersonic  phenomena  depending  on  body 
shape  dominate  the  flow  field  making  it  necessary  to  include  Mg^  in 
the  correlation,  and  the  characteristic  length  is  x^;  and  3)  Interpolation 
regime,  where  nondime nsionalization  of  the  time  dependent  Navier-Stokes 
equation  leads  to  the  normalizing  variables,  xtr^Ci^/Poo  and  Poo/PooT’ 
which  allow  a  single  correlation  of  the  interpolation  region  data. 

Only  in  a  vehicle  reentry  from  outside  the  atmosphere  down  to  sea  level 
can  one  take  a  single  vehicle,  make  a  single  test,  and  with  radar  follow  the  up¬ 
stream  motion  of  x{r  as  afunction  of  pressure  ratio  (altitude).  Suchinforma- 
tion  is  presently  unavailable  in  the  open  literature  and  so  no  comparison 
can  be  made  of  the  Map  predictions  with  an  appropriate  experiment 
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to  give  xtr  aa  a  function  ot  altitude  for  a  single  realization,  of  the  wake 
of  a  single  vehicle-.  The  comparison  of  the  Map  predictions  with  ballistic 
range  data  will  be  good  to  a  factor  of  two  because  all  ballistic  range 
data  available  to  date  was  used  to  derive  the  correlations,  good  to  a 
factor  of  tvv-o,  which  then  were  used  to  construct  the  Map  itself.  Hence, 
the  Map  certainly  applies  for  predicting  future  ballistic  range  results. 

In  ballistic  ranges,  each  firing  in  a  given  ambient  gas  produces  one  value 
of  x  and  this  value  will  be  predicted  by  the  Map  to  within  a  factor  of 
two. 

It  must  be  remembered  that  predictions  based  on  the  Transition 
Map  are  for  fluid  mechanical  unsteadiness  for  a  nonablating  body.  This 
fluid  mechanical  unsteadiness  is  observed  by  visual  analysis  of  schlieren 
and  shadowgraph  photographs.  Microwave  data  are  for  electron  wake 
unsteadiness.  It  is  most  likely  that  fluid  mechanical  unsteadiness  is  a 
prerequisite  for  electron  wake  unsteadiness.  Thus,  it  is  to  be  expected 
that  the  Map  predictions  are  closer  to  the  body  than  microwave  data  will 
indicate,  Fig; — Further,  the  selected  location  of  the  beginning  of  fluid 
mechanical  unsteadiness  as  determined  from  visual  examination  of  photo¬ 
graphs  is  not  sensitive  to  the  wave  length  of  the  unsteadiness  at  that 
point.  However,  the  location  of  wake  unsteadiness  as  determined  by 
microwave  scattering  is  very  sensitive  to  the  scale  of  the  electron  un¬ 
steadiness.  The  wake  after  having  gone  fluid  mechanically  unsteady  may 
not  produce  the  resonant  electron  fluctuation  scale  for  microwave  scat¬ 
tering  unti!  further  downstream  or  until  a  higher  ambient  pressure  is 
reached.  Thus,  at  a  given  pressure,  the  relative  location  of  the  radar 
wake  transition  further  downstream  than  the  fluid  mechanical  prediction 
is  a  reasonable  state  of  affairs,  even  without  considering  the  problem  of 
electron  production. 

The  work  here  for  the  hypersonic  wake  transition  map  has  been 
done  in  th*-  context  of  low  or  vanishingly  small  mass  transfer  from  the 
body  into  the  wake.  The  expectation  in  the  case  of  high  ablation  mass 
transfer  to  the  wake  would  be  the  following  for  each  respective  transi¬ 
tion  wake  regime:  1)  Far  Wake:  To  total -^G^A,  add  some  length  char¬ 
acterizing  the  rate  of  momentum  interchange  between  the  ablution 


.32 


/•(•fjritM'iw  hii'I  (mm  Mh-  riii  oimtu rr»d{  2)  Nn«i*  Wake;  Adjust  M^j  for  Jn» 

*  rpunwi  ll|irkiM  hh  at  Imumjary  layer,  lima  decreasing  and  mov¬ 

ing  th<?  t r *m nit  on  distance  forward,  ^ 


I'  lnrilly,  if*  the  precision  of  the  diagnostics  improve  and  ah  thw  depth 

of  understanding  increases  it  is  expected  that  modifications  in  detail  will 

be  necessary  for  the  results  presented  here.  For  example,  it  is  lively 

that  the  Near  Wake  transition  behavior  is  related  at  a  secondary  level  to 

the  width  of  the  wake.  Not  solely  dependent  on  distance  from  the  body  as 

43 

represented  in  the  Near  Wake  correlation  of  Eq.  (10).  Waidbuaser  and 
44 

Bailey  will  shortly  demonstrate  that  some  reduction  in  the  spread  of  the 
data  for  the  Near  Wake  correlation.  Fig.  13,  can  be  achieved  by  using 
•  x  •  6/d  as  the  ordinate,  where  6  is  the  wake  thickness,  instead 

of  just  p^  •  xfc  alone. 
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I  3  i ABSTRACT 

Hypersonic  Wake  Transition  data  over  a  range  of  Reynolds  numbers  have  been  obtained  over  the-past  several  years  at 
many  laboratories,  in  the  present  paper,  the  available  experimental  data  is  structured  into  three  laws;  The  Law  of  the  Far 
Wake,  The  Law  of  the  Near  Wake,  and  The  Law  of  die  Interpolation  Regime.  These  three  laws  are  then  represented  on  a 
single  Hypersonic  Wake  Transition  Map. 

Tin*  Reynolds  numb  r  is  selected  as  the  correlating  parameter,  where  local  conditions,  characterizing  the  state  of  the  gas 
at  transition  are  used  for  the  property  values.  The  correlating  constant  is  determined  from  the  ballistic  range  data.  For  hi 
Mach  number  flows,  binary  scaling  allows  the  representation  of  the  Reynolds  number  as  a  pressure  times  a  length. 

In  the  Far  Wake,  foi  the  first  onset  of  wake  unsteadiness,  the  wake  velocity  profiles  have  reached  their  asymptotic  behav¬ 
ior,  density  gradients  have  become  small,  and  details  of  body  shape  have  been  lost.  For  this  region,  the  characteristic 
length  is  represented  aW  O f>  A  based  on  the  total  drag  of  the  body.  The  Far  Wake  transition  rule  becomes  p  A/Cn  A 
—  constant.  rp  u 

In  thc.jsjcar  Wake,  It.cal  hypersonic  phenomena  depending  on  body  shape  dominate  the  flow  field.  Here  the  characteristic 
length  is  the  distant-,  from  ’.he  wake  origin  to  the  onset  of  unsteadiness,  xtr.  As  in  boundary  layer  flow  this  correlation 
shows  a  strong  dependence  on  an  appropriate  local  Mach  number,  M»,  at  the  point  of  wake  transition:  p  .  r,tr  *  constant, 
a  function  of  Me.  CO 

In  the  Interpolation  Regime,  one  traverses  from  the  Near  to  the  Far  Wake.  The  strong  favorable  pressure  gradient,  the 
strong  density  gradients  and  the  effect  of  body  shape  becomes  less  important  as  one  proceeds  downstream  along  the  Wake 
in  the  Interpolation  regime.  Thus,  ideas  derived  from  the  Navier- Stokes  equation  are  useful  and  lead  with  the  aid  of  the 
ballistic  range  data  t  a  sufficiently  accurate  general  interpolation  curve  to  connect  the  Ft  Wake. Law  with  the  Near  Wake 
Law.  The  Interpolation  curve,  like  the  p  .  VC p  A  Far  Wake  and  the  p  .  xtr  Near  Wake  correlation  is  good  to  a  factor 
of  2  in  the  pressure.  CO  00 

The  Transition  Map  provides  a  method  for  making  the  transition  curve,  xtr  as  function  of  pressure  ratio  (Altitude),  for  a 
given  body  at  a  giver,  flight  velocity:  Connect  the  Far  Wake  onset  altitude  for  theV  C  pA  (size)  of  the  body  with  the  Near 
Wake  ReXtr  behavior  for  the  corresponding  Msh  (shape)  of  the  oody.  Thus  a  complete  qualitative  Wake  Transition  Map 
is  presented  for  a  large  range  of  Mach  numbers  and  body  sizes  and  body  chapes. 

•This  paper  is  a  summary  of  the  author's  work  on  this  topic  over  the  several  vears  he  was  associated  with  A vco  Everett 
Research  Laboratory .  It  is  being  published  at  this  time  for  the  purpose  of  exchange  and  stimulation  of  ideas  and  not  as  an 
implicit  approval  of  the  findings  and  conclusions  therein. 
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